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Abstract 

Background: There is a strong need for early assessment of tumor response to chemotherapy in order to avoid the 
adverse effects of unnecessary chemotherapy and to allow early transition to second-line therapy. The purpose of 
this study was to determine the feasibility of ultrasonic spectral analysis for the in vivo characterization of changes 
in tumor microstructure in the evaluation of tumor response to chemotherapy using diagnostic ultrasound. 

Methods: Experiments were approved by the regional animal care committee. Twenty-four MCF-7 breast cancer 
bearing nude mice were treated with adriamycin or sterile saline administered by intraperitoneal injection. 
Ultrasonic radio-frequency (RF) data was collected using a clinically available ultrasound scanner (6-MHz linear 
transducer). Linear regression parameters (spectral slope and midband-fit) regarding the calibrated power spectra 
from the RF signals were tested to monitor tumor response to treatment. The section equivalent to the ultrasound 
imaging plane was stained with hematoxylin and eosin to allow for assessment of the density of tumor cell nuclei. 

Results: Treatment with adriamycin significantly reduced tumor growth in comparison with the control group 
{p = 0.003). Significant changes were observed in the ultrasonic parameters of the treated relative to the untreated 
tumors (p < 0.05). The spectral slope increased by 48.5%, from -10.66 ± 2.96 to -5.49 ± 2.69; the midband-fit 
increased by 12.8%, from -57.10 ± 7.68 to -49.81 ± 5.40. Treated tumors were associated with a significant decrease 
in the density of tumor cell nuclei as compared with control tumors (p < 0.001). 

Conclusions: Ultrasonic spectral analysis can detect changes in tumor microstructure after chemotherapy, and this 
will be helpful in the early evaluation tumor response to chemotherapy. 
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Background 

Tumor malignancy is one of the principal diseases that 
adversely affect human health and quality of life. More 
than half of all patients diagnosed with a malignant 
tumor will receive chemotherapy. At the present time, 
chemotherapy is still one of the most important cancer 
treatment methods. Early evaluation of tumor response 
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to chemotherapy in patients with cancer may help to 
avoid unnecessary treatment and enable the use of alter- 
native therapies. Currently in clinical oncology and ex- 
perimental therapeutics, assessment of tumor treatment 
response to chemotherapy relies on evaluating changes 
in tumor growth rate or volume weeks to months after 
the conclusion of a therapeutic protocol. These changes 
typically occur weeks to months late in the course of 
therapy. Functional techniques such as positron- 
emission tomography (PET), dynamic contrast-enhanced 
magnetic resonance imaging (MRI) and dynamic 
contrast-enhanced computed tomography (CT) have 
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been investigated regarding the early assessment of 
tumor response to chemotherapy; this involves depiction 
of reductions in the metabolic activity or the perfusion 
of the tumor, and some of the results have been promis- 
ing [1-3]. However, the use of such imaging modalities 
to monitor tumor responses to chemotherapy can be 
limited either by their cost or exposure of the patient to 
radiation. 

The aim of cancer therapy is to kill tumors by indu- 
cing cell death that can be used as an indicator of tumor 
response to therapy [4]. Increased tumor cell death early 
during the course of treatment, in preclinical and clinical 
studies, has been shown to be a good prognostic indica- 
tor of outcome [5,6]. Currently, standard methods for 
detecting cell death are invasive and require tissue bi- 
opsy for histologic analysis. Diffusion-weighted MRI has 
been used clinically to measure the increase in the ap- 
parent diffusion coefficient of water, which is thought to 
be increased in responding tumors owing to a decrease 
in cell density [7,8]. However, assessment of tumor re- 
sponse to chemotherapy requires repeated examinations 
and the cost of MRI limits its clinical use. 

Ultrasound is an attractive modality for the assessment 
of tumor response to therapy because of the ease with 
which it can be repeated without exposing the patient or 
animal to any risk of radiation. Ultrasound imaging sys- 
tems are also relatively inexpensive and mobile, a par- 
ticular benefit for animal studies. Ultrasonic spectrum 
analysis is used to extract information regarding tissue 
structures that is not conveyed in conventional B-mode 
imaging. Spectrum analysis of frequency-dependent 
backscattered radiofrequency (RF) data has been used to 
characterize tissue microstructures in the diagnosis of 
prostate cancer, ocular tumors, and cardiac abnormal- 
ities [9-11] and to differentiate benign lymph nodes from 
malignant lymph nodes [12]. It has been shown theoret- 
ically that the spectral parameters are related to tissue 
microstructural properties (e.g., effective acoustic scat- 
terer size and concentration) [13]. Ultrasonic spectral 
analysis with a high frequency transducer has been used 
to detect microstructural changes induced by radiother- 
apy, photodynamic therapy and chemotherapy in tumor 
xenografts [14-17]. Recently, Sadeghi-Naini and col- 
leagues used ultrasonic spectral analysis with a Sonix RP 
system at an approximate center frequency of 7 MHz to 
examine response to neoadjuvant chemotherapeutic in 
breast cancer patients and the preliminary results were 
promising [18]. However, this study was limited by the 
small number of cases (n = 10) and no confirmation of 
the changes in spectral parameters was obtained from 
evaluation of the corresponding tissue microstructures. 
The purpose of the present study was to determine the 
feasibility of ultrasonic spectral analysis for the in vivo 
characterization of tumor microstructural changes in the 



evaluation of tumor response to chemotherapy using 
diagnostic ultrasound. 

Methods 

Animal model 

This study was approved by the Committee on the 
Ethics of Animal Experiments of the Sun Yat-Sen 
University under the guidelines of the National Institutes 
of Health for the care of laboratory animals. Human 
breast cancer cell line MCF-7 was obtained from State 
Key Laboratory of Oncology in Southern China. MCF-7 
cells were grown in DMEN culture medium (Hyclone 
Co., UT, USA) supplemented with 10% fetal bovine serum 
(Gibco, Grand Island, NY, USA), penicillin (50 U/ml), and 
streptomycin (50 (ig/ml) at 37°C in a humidified 5% C0 2 
atmosphere. For inoculation, approximately 4 x 10 8 MCF- 
7 cells suspended in phosphate-buffered saline were 
injected subcutaneously into the right axillary fossa of 24 
5- to 6-week old BALB/c nude female mice. 

Twenty-four nude mice were randomly divided into 
two groups. Twelve of the nude mice were treated with 
adriamycin (Shenzhen Main Luck Pharmaceuticals Inc., 
Guangdong, China) diluted in sterile saline and adminis- 
tered once daily by intraperitoneal injection (4 mg/kg). 
Drug administration began at day 10 post tumor cell im- 
plantation when tumors had reached a size of 5 mm. 
The remaining twelve nude mice in the control group 
were given the vehicle control medium (sterile saline) 
according the timing and dosing schedule used for the 
treated group. 

Ultrasound data acquisition 

Ultrasound imaging was performed at day 18 (i.e. 7 days 
after initiation of therapy). For the ultrasound imaging 
studies, each mouse was anesthetized by intraperitoneal 
injection of pentobarbital sodium at a dose of 75 mg/kg 
(Sigma, St. Louis, MO, USA). Centrifuged gel was used 
to minimize bubble formation in the gel and a stand-off 
gel pad was placed on the skin for scanning. A com- 
mercially available clinical ultrasound scanner, Sonix 
TOUCH (Ultrasonix Medical Corporation, Richmond, 
Canada) with a 6-MHz linear transducer was used to 
simultaneously collect B-mode images and RF data 
from the treated and control tumors. For data acquisi- 
tion, the ultrasound transducer was positioned such 
that the focal zone was at the same depth in each im- 
aged specimen to control for any potential attenuation. 
All RF data were sampled with 16 bit resolution at a 
frequency of 35 MHz. All images and radiofrequency 
data were digitally recorded. All ultrasound examina- 
tions were performed by one radiologist, who was 
blinded to the treatment status. The greatest longitu- 
dinal, transverse and anteroposterior dimensions of tu- 
mors were measured in fundamental grayscale imaging 
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using calipers. Tumor volume was calculated using the 
formula for a prolate ellipsoid: volume = tt/6 x length x 
width x depth. The largest cross-section plane of the 
tumor was imaged with the transducer held manually 
in this position throughout the examination. To evalu- 
ate the echogenicity changes in the tumor after chemo- 
therapy, Adobe Photoshop 6.0 (Adobe Systems, San 
Jose, CA, USA) was used to measure the gray scale in- 
tensity of the ultrasound images. 

Ultrasonic spectral analysis 

The ultrasound RF data from was imported into 
MATLAB-based (v. 2009a: MathWorks, Natick, MA, 
USA) software developed in our lab for ultrasound spec- 
tral analysis. For the tumor, rectangular regions of inter- 
est (ROI) were centered approximately at the focal 
depth of the transducer. Three representative ROIs were 
selected for each tumor sample and averaged for the 
final analysis. RF data from each line segment were 
multiplied by a Hamming weighting function to sup- 
press spectral lobes and the Fourier transform was 
computed. 

The power spectrum was obtained by averaging the re- 
sults from the independent scan lines. This power 
spectrum was divided by the power spectrum of the 
echo from a calibration target. A quartz flat was used as 
the calibration target and the perpendicular reflection 
off the quartz flat located at the focal point of the trans- 
ducer was used to derive the power spectrum. Linear re- 
gression analysis was applied to the calibrated spectral 
amplitude to provide a best-fit line (Figure 1) [19]. To 
assess the role of ultrasonic spectrum analysis in moni- 
toring tumor response to chemotherapy, two spectral 
parameters from the regression analysis were computed: 
the spectral slope (SS) and the midband fit (MBF). The 
SS is the slope of the linear regression of the calibrated 
spectrogram and the MBF is the value of the regression 
fit at the center frequency over which the spectrum was 
measured. The spectral parameters can potentially be 
mathematically related to the physical characteristics of 



ultrasound scatterers, such as the size and concentration 
of the ultrasound scatterers [20]. Further details on the 
theoretical and signal analysis considerations and the re- 
lationship between spectral parameters and tissue micro- 
structure can be found elsewhere [19,21]. 

Histopathologic examination 

At the end of the experiment, mice were sacrificed using 
the standard method. Tumors were removed and fixed 
in 10% buffered formalin before paraffin processing. The 
tumor specimens were sectioned (5 \im) at the largest 
cross sections corresponding to the ultrasound imaging 
planes. Sections were stained with hemotoxylin and 
eosin (H&E) and assessed microscopically for changes in 
cell morphology. 

Regions with the highest density of tumor cell nuclei 
were located by scanning the tissue sections under a 40x 
power microscope. After identification of the regions of 
highest density, ten different fields were randomly 
chosen within these regions at 400x power. The 400x 
histology images were then analyzed to measure the 
density of tumor cell nuclei by counting the number of 
nuclei in each image using Image Pro Plus software 
(Image Pro Plus 6.0, Media Cybernetics, Silver Spring, 
MD, USA) (Figure 2). The H&E stained image segmenta- 
tion was based on the following Hue-Saturation-Intensity 
(HSI) parameters: Hue (0-255), Saturation (0-255) and 
Intensity (0-120). The segmented areas in the images 
were filtered to count blue nuclei. This filtering used 
thresholds as follows: area (minimum = 50 pixels) and box 
x/y (minimum = 0.5; maximum =2). The "split objects" 
function was used to separate cells touching each other. 
The average count of ten 400x images was used for the 
statistical analysis. 

Statistical analysis 

All statistical analyses were performed using SPSS version 
16.0 (SPSS, Inc, Chicago, IL, USA). The Kolmogorov- 
Smirnov test was applied for evaluation of normal distri- 
bution and the Levene test for evaluation of homogeneity 




B 



O 
3 



s 

< 



•50 



-150 



Calibrated Spectrum 
Linear Regression 




Frequency (MHz) 



Figure 1 Ultrasound image and power spectrum. (A) Ultrasound image of the MCF7 tumor with the ROI located in the subcutaneous tumor 
tissue and (B) the corresponding calibrated power spectrum (blue curve) with linear regression (red line). 
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Figure 2 Histologic section and binary image. (A) A hemotoxylin and eosin (H&E) stained image of the MCF-7 tumor and (B) the 
corresponding binary image indicating the presence of nuclei were identified using Image Pro Plus. The red square showed the area threshold 
used for counting the nuclei. 



of variance. Student s £-test was used to determine the sig- 
nificant differences in measurement data between the 
treated and control tumors. A leave-one-out cross- 
validation method was used to test the ability of slope and 
midband-fit to distinguish treated tumors from control tu- 
mors. A p value of 0.05 or less was considered as being 
statistically significant. 



Results 

Effect of adriamycin on tumor growth 

After treatment for 7 days, the mean tumor volumes of 
the treated and control tumors were 0.08 ± 0.03 cm 3 and 
0.17 ± 0.08 cm 3 , respectively. Treatment with adriamycin 
(4 mg/kg once daily) significantly reduced tumor growth 
in comparison to the control group (p =0.003). 



Changes in ultrasonic spectral parameters 

After treatment for 7 days, the mean gray scale intensity of 
the control and the treated tumors in conventional B-mode 
ultrasound images was 15.40 ± 3.86 units and 19.00 ± 1.25 
units, respectively. Treatment with adriamycin (4 mg/kg 
once daily) significantly increased the gray scale intensity of 
the conventional B-mode ultrasound images as compared 
with control tumors (p = 0.009). Spectrum analysis of 
frequency-dependent backscattered radiofrequency data 
showed the difference in spectral parameters between the 
treated and control tumors. The mean spectral slope of the 
control and treated tumors was -10.66 ± 2.96 dB/MHz 
and -5.49 ± 2.69 dB/MHz, respectively. Seven days 
treatment with adriamycin (4 mg/kg once daily) signifi- 
cantly increased spectral slope (by 48.5%) in comparison to 
the control group (p < 0.001). The mean MBF of the 
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Figure 3 Representative ultrasound images and corresponding ultrasound spectral parameters. Ultrasound images for (A) a control tumor 
and (B) a treated tumor, and the corresponding ultrasound spectral parameter (C) characterization of tumor responses to chemotherapy. The 
ultrasound images showed a noticeably brighter grayscale intensity in treated tumors (B) relative to control tumors (A). (C) Ultrasonic spectral 
analysis indicated a separation of the regression line between the control and the treated tumors. 
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control and treated tumors was -57.10 ± 7.68 and -49.81 ± 
5.40 dB, respectively. MBF was significantly increased (by 
12.8%) after 7 days when the treatment tumor group was 
compared with the control group (p = 0.013) (Figure 3). In 
distinguishing between treated and control tumors using 
"leave-one-out" cross-validation, the correctly classified 
rate of slope and MBF was 87.5%, and the sensitivity and 
specificity were 83.3% and 91.7%, respectively. 

Histological changes 

The most prominent histological changes after chemo- 
therapy in the treated tumors were related to the density 
of tumor cell nuclei, nuclear size and the extent of cyto- 
plasmic and nuclear vacuolation of the tumor cells. The 
number of tumor cell nuclei evaluated in histological 
slice from treated tumors was 78.51 ±13.11 counts per 
high-power field (HPF) and in control tumors was 
334.50 ± 44.57 counts per HPF. Treatment with adriamycin 
(4 mg/kg once daily) significantly reduced the density 
of tumor cell nuclei in comparison to the control group 
(p < 0.001). H&E staining revealed other microstruc- 
tural change in treated tumors, involving nuclear struc- 
ture manipulation (condensation and fragmentation) 
(Figures 4 and 5). 



Correlation between ultrasonic spectral parameters and 
histological results 

The density of cell nuclei density was found to be nega- 
tively correlated with spectral slope (r = -0.670, p < 0.001) 
and MBF (r = -0.450, p = 0.027). 

Discussion 

Reduction in tumor size is a late sign of effective chemo- 
therapy, early response to chemotherapy is difficult to 
assess using conventional radiographic modalities [22]. 
Currently many cancer treatment regimens require sev- 
eral courses of chemotherapy before it can be deter- 
mined whether or not the treatment has been effective. 
The availability of non-invasive methods for predicting 
and/or detecting therapeutic response to chemotherapy 
at an early stage of treatment would facilitate the ra- 
tional design and individualization of therapy protocols 
for cancer patients and allow early transition to second- 
line therapy. The present study was undertaken to assess 
the usefulness of ultrasonic spectral analysis in the 
evaluation of conventional chemotherapy in a murine 
breast cancer model using tumor volume and the dens- 
ity of tumor cell nuclei as gold standards. This study has 
shown for the first time, using a clinically available 6 





Figure 4 Histopathologic analysis of tumor cell nuclei density. (A) The graph displays the mean density of tumor cell nuclei in the control 

and treated tumors (* = p < .001). (B) and (C) Representative photomicrographs of hemotoxylin and eosin stained sections of the control (B) and 

treated (C) tumors. Staining revealed microstructural changes in the tumor treated with adriamycin, including a decreased density of cell nuclei, 

and cytoplasmic and nuclear vacuolation, and clumping of nuclear chromatin (original magnification, x400). 
\ J 
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Figure 5 Scatter plots of changes in spectral parameters. (A) Spectral slope, (B) midband-fit versus cell nuclei density and (C) spectral slope 
versus midband-fit. 



MHz ultrasound transducer, that chemotherapy effects 
can be characterized by means of ultrasonic spectral 
analysis in preclinical mouse cancer models. The 
changes in spectral parameters were interpreted as oc- 
curring as a direct consequence of cell death and 
changes in the density of cell nuclei after chemotherapy. 
This method was able to detect changes in the solid 
tumor microstructure after chemotherapy and, conse- 
quently has the potential to detect tumors that are re- 
sponsive to treatment earlier than using conventional 
methods. 

In the current study, histologic analysis revealed that 
the most prominent microstructural changes after 
chemotherapy were related to the density of cell nuclei, 
cytoplasmic and nuclear vacuolation and clumping of 
nuclear chromatin as indicated by the H&E staining; 
these changes are frequently seen in breast cancer fol- 
lowing neoadjuvant chemotherapy [23,24]. Our histo- 
logic observations suggested that the main changes in 
the ultrasonic spectral parameters after chemotherapy 
were related to microstructural changes regarding both 
the density of cell nuclei and cytomorphologic changes. 

It has long been suggested that ultrasound backscatter 
variables may relate to the physical properties of tissues 
[20]. Ultrasonic spectral analysis techniques have been 
used by many investigators to add information regarding 
tissue microstructure to images generated by conventional 



ultrasound [25]. The spectral slope is an indicator of ef- 
fective scatterer shape and size and an increase in spectral 
slope corresponds to a decrease in effective scatterer size 
[26]. The observed changes in the spectral slope after 
chemotherapy could be suggestive of structural changes in 
the tumor cells after treatment, namely cell shrinkage, nu- 
clear condensation and fragmentation. The MBF is an- 
other measure of ultrasound backscatter and depends on 
multiple factors, including scatterer shape, size, concentra- 
tion, and the acoustic impedance change between the 
acoustic scatterer and the surrounding medium [26]. In- 
crease in MBF was seen in treated tumors after chemo- 
therapy relative to control tumors. The mechanism 
behind this increase was broadly linked to changes in cell 
and nuclear morphology observed histologically after cell 
death. The strength of ultrasound backscatter depends 
both on scatterer size and differences in the mechanical 
properties of the scatterer, the surroundings (compress- 
ibility and density) and the scatterer number density (e.g., 
how many scatterers there are per unit volume) [15]. The 
changes in cell and nuclear sizes during the sequence of 
cell death resulted in an increase in the variance of cell 
sizes; a previous study by Vlad et al. [27] demonstrated 
that an increase in cellular size variance contributes to the 
increase in ultrasound backscatter during cell death. The 
decrease in the density of tumor cell nuclei and extensive 
changes in cytoplasmic and nuclear vacuolation of the 
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tumor cells might result in changes in acoustic impedance, 
contributing to the increase in MBF. 

Ultrasonic parameters, including spectral slope and 
MBF, have previously been used to characterize diseased 
tissue or tissue and cell samples exposed to different 
therapeutic agents [14-17,27,28]. Kolios and colleagues 
used ultrasonic spectral analysis to measure changes in 
the spectral slope and MBF in vitro for cell samples ex- 
posed to chemotherapeutic drugs [29]. Significant in- 
creases in spectral slope and MBF in treated cells were 
observed after exposure to chemotherapy that were in 
close agreement with theoretical predictions. However, 
more complex histological changes were observed in 
in vivo in mouse breast cancer models after chemother- 
apy than in vitro in cell samples exposed to chemothera- 
peutic drugs. In our study, other than cell shrinkage and 
nuclear condensation and fragmentation, the most prom- 
inent microstructural changes observed after chemother- 
apy were related to the decrease in the density of tumor 
cell nuclei, and cytoplasmic and nuclear vacuolation as re- 
vealed using H&E staining. Moreover, high-frequency 
ultrasound has been used in previous studies to detect 
apoptosis in vivo in cell samples and animal systems ex- 
posed to different anticancer therapies [14-17,29]. A pene- 
tration depth of 2-5 cm at frequencies of 10-30 MHz 
limits the applicability of this technique to superficial 
regions. The present study confirmed that ultrasonic 
spectral analysis using a lower frequency of 6 MHz 
could be used to monitor changes in tumor microstruc- 
ture in a mouse breast cancer tumor model after 
chemotherapy. Sadeghi-Naini and colleagues used 
ultrasonic spectral analysis with a Sonix RP system at 
an approximate center frequency of 7 MHz to examine 
neoadjuvant chemotherapeutic response in breast can- 
cer patients, and their preliminary results were promis- 
ing [18]. Lower frequency ultrasound penetrates much 
deeper, and previous studies have confirmed that the 
backscattered data from low frequency ultrasound 
could be potentially subjected to the same analysis to 
provide information regarding structural changes in the 
tissue at the cellular level [30-33]. 

Increased tumor cell death occurring early during the 
course of treatment, in both preclinical and clinical stud- 
ies, has been shown to be a good prognostic indicator of 
outcome [5,6]. As a result of tumor cell death, the dens- 
ity of tumor cell nuclei will decrease significantly after 
chemotherapy [23,24,34]; this has become a marker for 
the detection of early indications of tumor response to 
chemotherapy [4], Noninvasive diffusion-weighted MRI is 
a well studied imaging technique for quantifying the in- 
crease in the apparent diffusion coefficient of water caused 
by a decrease in tumor cell density within 2-4 days, prior 
to visible changes in tumor morphology or size in patients 
with breast [35], brain [36], and ovarian [8] cancers that 



responded to treatment. However, the cost of diffusion- 
weighted MRI limits its extensive use in monitoring tumor 
response to cancer therapies. In contrast, in the measure- 
ment of changes in water diffusion caused by cell death, 
ultrasonic spectral analysis provides a simple imaging ap- 
proach that can directly characterize tumor microstruc- 
tural changes (decreased density of tumor cell nuclei) 
after chemotherapy. 

A potential limitation of our technique was that there 
were potential variations in matching the ultrasound 
image planes with the histological slices, due to differ- 
ences in slice thickness between the ultrasound scanner 
and histological sections. In order to make ultrasound 
data correspond more closely to histologic measure- 
ment, the largest cross-section planes were used in both 
techniques. Another limitation was that there was no 
compensation for attenuation in our study, However, be- 
cause of the fact that the tumors were close to the skin 
surface, and the center of the ROI was located at ap- 
proximately 0.2 cm below the skin surface, attenuation 
compensation will not likely lead to significant changes 
in the results for a 6 MHz transducer. Moreover, the 
limiting scatterer size for a 6-MHz linear transducer will 
be about 85 \xm (13), which is much larger than the nu- 
clear diameter (about 15 \im); consequently, it was diffi- 
cult to identify the histological texture of the scatterer 
for the 6-MHz transducer, and therefore to specify the 
histological changes that caused the changes in the spec- 
tral parameters. 

Conclusions 

In conclusion, this study indicated that ultrasonic spec- 
tral analysis provided a simple way to characterize tumor 
microstructural changes after chemotherapy using a 
clinically available 6 MHz ultrasound transducer. A sig- 
nificant increase in spectral slope and MBF were 
detected using this noninvasive imaging technique after 
chemotherapy. This would allow tumor imaging before 
and at multiple times during treatment without the need 
for injecting specialized contrast agents as is required 
using other techniques (e.g. PET, dynamic contrast- 
enhanced CT and dynamic contrast- enhanced MRI). This 
noninvasive technique shows considerable potential in the 
early assessment of tumor response to chemotherapy. 
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